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Q^ '■ Abstract 

p ,' Corrections ~ a^ to the positronium decay rate, induced by one-loop 

^, ■ vacuum polarization diagram, are calculated. Their relative values are 



v2 r 1 r^ r^r.n/o\/ / \2 



0.4468(3) (a/7r) for para- and 0.960(3) (a/7r) for orthopositronium. 



The measured value of orthopositronium decay rate is Hj 

Te^p = 7.0482(16)/is-\ 

At the same time, its theoretical value (including corrections ~ (a/7r) and 
~ a^ log a) is firmly established [@, |^, g, |5], g, [7| and constitutes 



ri,k = ma" ?<^^ [l - 10.28^ - ^a^ log - 

2 



7.03830/is-\ 



To lift the disagreement between them by (a/Ti) -corrections, the factor 
at {a /it) should be 250(40), which may look unreasonably large and very 
unusual for QED. To remove this discrepancy would be of importance for 
clear understanding of the relativistic bound-state problem in QED. 

However, there exists at least two contributions of the second order with a 
large coefficient at {a/ir) . The first one arises when we square the one-loop 
correction to the annihilation amplitude [^ , it is equal to 

r2th = 28.8(2)(a/7r)2ro = 0.00112(l)/is-^ 

(here Fq is the lowest-order decay rate ). The second one is the relativistic 
correction calculated in 



r^th = 46(a/7r)2ro = 0.00179(12)/is"^ 

Calculation of the contributions of two-loop corrections to the annihilation 
amplitude is required to obtain the total correction ~ a^ to the orthopositro- 
nium decay rate. There exists an opinion [|TDj that numerical factor at (a/7r) 
in the amplitude can also be large (actually, the above-mentioned relativistic 
correction 191 is contained also in that contribution). 



There is a large number of diagrams which contribute to the second- 
order amplitude^, but only sum of them is gauge invariant. Nevertheless, 
Tf at the calculation the gauge does not contained infrared singularities 



there is a class of the two-loop diagrams sum of which is gauge invariant 
by itself. It consist of diagrams which contain photon self-energy insertion 
( these diagrams for orthopositronium (o-Ps) and parapositronium (p-Ps) 
are displayed in Fig. 2 and 1 respectively). In this paper we calculate these 
corrections both for o-Ps and p-Ps. 

To account for the self-energy part in photon propagators, we substitute 

' ' + 1 Ri^) i^r^—d^ . (1) 



(Kallen — Lehmann representation); to first order in a, 



a 



s — Am?' s + 2m^ 



Hence the considered contribution to positronium decay rate is of the form 

^^ = IL ^(^) ^(^) ^^ ' 

where Q{s) is the one-loop correction to the decay rate calculated with the 
massive virtual photon (m^ = ^/s ). 

The calculation of corrections ~ a^ that are induced by vacuum polariza- 
tion, was performed in the same manner as calculation of corrections ~ a, 
which is described in the previous papers 0, 0, §| . Our calculation dif- 
fers from the one-loop one by the substitution ([I|) , in which we must deal 
only with second term in the right-hand side, that introduced a one more 
integration compared to one-loop case. Appearing integrals were found 
numerically using convenient substitution s = 4m^/(l — x^) . To sufficient 
accuracy the initial particles may be considered being at rest. We use co- 
variant summation over polarizations of the final photons. The results of 

is used, there are 102 diagrams for orthopositronium and 36 diagrams for 
parapositronium. 



computing of various diagrams are presented in Tables 1,2. Thus, the total 
corrections to the positronium decay rate, induced by vacuum polarization, 
constitute 

^r?/' = 0.4468(3) {a/irf T^^' = 0.0194^8-^ , 

5Tl-^' = 0.960(3) {a/nf T'l^^' = 3.735(12) • lO^V"^ • 

Let us note that the main part of the complete result originates as con- 
tribution of the diagram 2(f), which equals to — |(a/7r)r^ ( here F^ is part 
of a-correction, which arises due to graph shown in Fig. 3). 

In order to check the results of the numerical calculation we also com- 
puted differential decay rate of o-Ps using three-dimensionally transverse 
summation over the final photons polarizations. The obtained result coin- 
cides with the differential width which was calculated employing the covari- 
ant photon polarization sums ( though contributions of individual graphs do 
depend on manner of the summation, with the exception of the contribution 
of 2(f), which is gauge invariant by itself). 

Another method to check our results is as follows. Let Ti{s) to be the 
contribution to the Q{s) due to the diagram [i) [i = a,b,c for p-Ps and 
i = a, ..,f for o-Ps ). Their asimptotics at s ;> 1 are Ti{s) = Ai/s. The 
values of the Ai^s were calculated analytically for all contributions to the 
p-Ps width and for i = a,b,c,d in the case of o-Ps. We examined that the 
numerically obtained T^'s have the asimptotics which coincide with those 
found analytically. In the case of o-Ps the check was made for the differential 
width. In addition, for the graph 2(/) Af may be related to the part of 
correction ~ a to the o-Ps decay rate, which arises due to graph 3; Af 
which was found in this way is in agreement with the numerical results. 

If we add obtained in that paper contribution to the orthopositronium 



decay rate FlJi^ * to ones calculated earlier Fi^/j, T2th', ^sth^, it is obtained 
the result which still too different with the experimental one. It is conceiv- 
able that further two loops calculations will result in the better accordance 
between theory and experiment. 

One of the authors (A.B.) wish to thank I.B.Khriplovich for useful dis- 
cussions. 



Table 1. Various contributions to p-Ps decay rate. 



a 


h 


c 


total 


0.0158(1) 


0.0975(1) 


0.3335(1) 


0.4468(3) 



Table 2. Various contributions to o-Ps decay rate. 



a 


h 


c 


d 


e 


/ 


total 


0.151(1) 


-0.0895(3) 


0.210(1) 


0.062(1) 


-0.094(1) 


0.720(1) 


0.960(3) 
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Fig. 1 
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Fig. 2 
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